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Abstract 



A partial-wave analysis of all pp scattering data below 925 MeV/c antiproton 
laboratory momentum is presented. The method used is adapted from the 
Nijmegen phase-shift analyses of pp and np scattering data. We solve the 
Schrodinger equation for the coupled pp and ran channels where the long- and 
intermediate-range interactions are described by a theoretically well-founded 
potential. This gives the rapid variations of the scattering amplitudes with en- 
ergy. This potential consists of the Coulomb potential with the main relativis- 
tic correction, the magnetic-moment interaction, the one-pion-exchange po- 
tential, and the heavy-boson exchanges of the Nijmegen one-boson-exchange 
potential. Slow variations of the amplitudes due to short-range interactions, 
including the coupling to mesonic annihilation channels, are parametrized by 
an energy-dependent, complex boundary condition, specified at a radius of 
r = 1.3 fm. The Nijmegen 1993 pp database, consisting of 3646 pp scatter- 
ing data, is presented and discussed. The best fit to this database results 
in Xmin/-^data = 1-043. This good fit to the data shows that the Nijmegen 
long- and intermediate-range potential is essentially correct. The pseudovec- 
tor coupling constant of the charged pion to nucleons is determined to be f£ 
= 0.0732(11) at the pion pole, where the error is statistical. 
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I. INTRODUCTION 



A partial-wave analysis (PWA) of all antiproton-proton (pp) scattering data below an- 
tiproton laboratory momentum p\ a ^ = 925 MeV/c is presented. In the field of nucleon- 
nucleon (AW) scattering, phase-shift analyses (more properly called partial-wave analyses, 
PWAs) have a long history and at present the multienergy or energy-dependent partial- wave 
analyses of AW scattering data have reached a stage of considerable sophistication [0-0]. Due 
to the poor quality of low-energy antiproton beams and the resulting absence of accurate ex- 
perimental data, analogous model-independent studies of the much more complex pp system 
have in the past always been impossible. 

In recent years, however, experimental progress has been very significant, in particular 
due to the advent in 1983 of the LEAR facility at CERN. In the momentum region that we 
consider in the analysis presented in this paper, the situation between 400 and 925 MeV/c is 
quite good: a variety of observables have been measured with impressive accuracy. However, 
the practical difficulties involved in constructing a high-quality antiproton beam of even lower 
momentum are large. As a result, the pp database below about 400 MeV/c is still by far 



not as good as one would like PHlDl, in contrast to the case of AW scattering, where for 
instance very accurate proton-proton differential cross sections have been taken at energies 
as low as T lab = 0.35 MeV (p lah = 25.6 MeV/c). 

During the last 10 years a new method has been developed by the Nijmegen group to 
perform partial-wave analyses of the AW (proton-proton and neutron-proton) scattering 
data below laboratory kinetic energy = 350 MeV (p\ a b — 883 MeV/c). The hallmark 
of this method is that the theoretical knowledge of the AW interaction is exploited as much 
as possible in the description of the energy dependence of the partial-wave scattering am- 
plitudes. This is done by solving the relativistic Schrodinger equation for each partial wave 
and at each energy with the theoretically well-known long-range AW interaction represented 
by the potential for r > b = 1.4 fm. This long-range interaction is responsible for the 
fast variation with energy of the scattering amplitudes. The much slower energy variations 
of the amplitudes due to the short-range interactions are parametrized phenomenologically 
by the energy dependence of a boundary condition at r = b. In this way one also avoids 
the complications due to the lack of knowledge of this short-range interaction. In order 
to achieve a good fit to the AW data one has to include in the long-range interaction the 
complete electromagnetic interaction (including relativistic and some two-photon-exchange 
effects, the magnetic-moment interactions, and the vacuum-polarization potential), the one- 
pion-exchange potential, and also the intermediate-range forces (due to two-pion exchange 
and/or heavier-meson exchange) of some realistic AW potential model. 

The dominant feature of pp scattering at low energies is the annihilation into mesons, 
a process that has no counterpart in AW scattering. Annihilation is, in principle, a very 
complicated spin-, isospin-, and energy-dependent multiparticle process. At rest, where 
kinematically the production of 13 neutral pions is allowed, 5 pions are produced on the 
average |TT|], and of the order of 100 two-meson channels contribute significantly |T^,T3]. In 



recent years some progress has been made in understanding specific annihilation processes 
in terms of quark-gluon degrees of freedom. However, for the description of elastic pp —>■ pp 
and charge-exchange pp —>■ nn scattering, only a phenomenological approach to annihila- 
tion is feasible at present. In potential models for antinucleon-nucleon (AW) scattering one 
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usually simplifies things drastically by taking the annihilation potential to be completely 
independent of spin, angular momentum, isospin, and energy. This assumption is imple- 
mented either by applying a simple absorptive boundary condition fI3}-(r7||, or by using a 
state-independent two- or three-parameter optical potential [p^ -p3|. However, when one is 
interested in describing the data quantitatively, including spin-dependent observables, a less 
naive approach is called for. For instance, in the Paris AW model p4| , p5| a spin-, isospin-, 
and energy-dependent optical potential is employed, and in the Nijmegen coupled-channels 
model p6[ each AW channel is coupled to two effective two-body mesonic channels (for the 
coupled-channels approach see also Ref. J27]]). In both cases of the order of 15 parameters 
were needed to obtain a more-or-less satisfactory fit to the pre-LEAR data, the bulk of 
which consisted of elastic differential cross sections. At that time charge-exchange data and 
spin-dependent observables were practically absent. Because the NN system is so much 
more complicated than the already quite complex NN system it was then believed that it 
would be impossible to perform a PWA of the pp scattering data. 

The complexity of the pp system when compared to the NN system (below the pion- 
production threshold) is reflected in the following manner in a PWA. In proton-proton (pp) 
scattering (isospin 1 = 1) one has to specify at each energy 2 phase shifts ( 1 5o and 3 Po) 
for J = and on the average 2.5 phase parameters for each value of the total angular 
momentum J ^ 0. As an example: for J = 1 one needs only one phase shift ( 3 P\), while for 
J = 2 one needs 4 phase parameters (the 1 D 2 , 3 P 2 , and 3 F 2 phase shifts and the e 2 mixing 
parameter). In a PWA of the neutron-proton (np) scattering data (both isospin 1 = and 
1 = 1) again 2 phase shifts are required for J = 0, but now 5 phase parameters are required 
for each value of J ^ 0. Due to the lack of sufficient high-quality np data it has been 
impossible to do a good PWA of the np data alone. One needs to take the 1 = 1 phases 
(with the exception of the 1 Sq phase shift) from the PWA of the pp data and correct them 
for electromagnetic and mass- difference effects (M p versus M n and m^o versus m^±). In the 
case of pp scattering the Pauli principle is not operative. Apart from this, the possibility for 
the pp system to annihilate into mesons complicates things even further. This means that 
one has to determine no less than four times as many phase parameters compared to the 
case of np scattering, so 8 phase parameters are required for J = and 20 phase parameters 
for each value of J ^ 0. In view of this, the situation with regards to a pp PWA indeed 
seemed quite hopeless. 

Using essentially the same strategy as in the Nijmegen multienergy partial-wave analyses 
of AW scattering data, and with the available recent high-quality data from LEAR and KEK, 
we have nevertheless been able to perform a PWA of pp scattering data below 925 MeV/c. 
This work was started in 1987 and a preliminary report of this analysis has already been 
given [p9[| . The Schrodinger equation for the coupled antiproton-proton and antineutron- 
neutron (nn) channels is solved. The short-range interaction, including the coupling to 
the mesonic annihilation channels, is parametrized by way of a complex boundary condition 
specified at r = b = 1.3 fm. The long-range interaction consists of the Coulomb potential, the 
magnetic-moment interaction, and the one-pion-exchange potential. The tail of the heavy- 
boson-exchange part of the Nijmegen potential |30] is used as intermediate-range interaction. 
A lot of time and effort has gone into collecting, scrutinizing, and cleaning up the world set 
of pp scattering data, which contains quite some flaws and contradictory data. Exactly the 
same arguments were used in this process as were used in the set-up of the Nijmegen AW 



3 



database |5|-[7|]. The resulting Nijmegen 1993 pp database in the momentum interval 119-923 
MeV/c consists of iV data = 3646 pp data, which can be fitted with Xmm/^data = 1.043. In 
view of the excellent quality of this fit one can be confident that most of the amplitudes are 
quite well described. The same methods have also been applied by us to the strangeness- 
exchange reaction pp —>■ AA, which is an even more complicated process pi,52|. 

This paper is organized as follows. In Sec. II we present the details of the method of 
analysis. In Sec. Ill the treatment of the short-range interaction by parametrizing an energy- 
dependent boundary condition is discussed. In Sec. IV the long-range electromagnetic and 
pion-exchange potential are presented and we specify the heavy-meson exchanges used as 
intermediate-range interaction. In Sec. V we discuss the non-trivial problem how to extract 
the nuclear scattering amplitude in the presence of electromagnetic effects. The definition 
of the phase-shift and mixing parameters for antinucleon-nucleon scattering can be found 
in Sec. VI. The statistical tools used in the analysis are briefly reviewed in Sec. VII. In 
Sec. VIII the Nijmegen 1993 antiproton-proton database is extensively discussed. Next, 
in Sec. IX we present the results of the analysis, including the determination of the pion- 
nucleon coupling constant. The most important results and conclusions are summarized in 
Sec. X. The algorithm to extract the phase parameters from the S matrix is reviewed in the 
Appendix. 



II. THE METHOD OF ANALYSIS 



The two-body scattering process is described with the coupled-channels relativistic 
Schrodinger equation 



[A +p 2 - 2mV] ip(r) = 



(1) 



This is a matrix equation in channel space. We use the physical particle basis pp, nn in order 
to treat electromagnetic effects properly and to account for the threshold of charge-exchange 
scattering pp — > nn at p^y, = 99.1 MeV/c (Ti ab = 5.2 MeV). The connection between the 
channel momentum p and the total energy y/s in the center-of-mass system is given by 



the relativistic expression p 



is 



m 



(for equal masses). The relativistic Schrodinger 



equation Eqn. (JJ) is a differential form of the relativistic Lippmann-Schwinger integral 
equation [|33,34"|. The difference between the relativistic and the ordinary non- relativistic 



Lippmann-Schwinger equation |35] is the relation used between energy and momentum. The 
relativistic Lippmann-Schwinger equation is in turn equivalent to three-dimensional relativis- 
tic integral equations like the Blankenbecler-Sugar equation |36|-p9^ . For a discussion about 
the derivation of potentials for use in the relativistic Schrodinger equation, starting from 



see for instance Refs. 42.13 



the field-theoretical Bethe-Salpeter equation [f40| , [4T 

The interaction in the region r > b is described by a theoretically well-founded 
antinucleon-nucleon potential. This potential is given by 



V = V c + Vmm + Vi 



N 



(2) 



where Vc and Vmm are the relativistic Coulomb and magnetic-moment interaction respec- 
tively, is the NN meson-exchange potential. The precise forms of these potentials are 
discussed in Section IV. After making the partial-wave projection by writing 
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we obtain the radial Schodinger equation, which for partial waves with total angular mo- 
mentum J reads 



dr 2 



+ P 



2mV- J 



. 



(5) 



V J is now a matrix with elements (£'s'a'\V J (r)\£ s a), where a is an index used to distinguish 
the different channels. For partial waves with i = J, s = 0, 1 all matrices are 2x2, and for 
waves with £ = J ±1, coupled by the tensor force, they are 4x4. The Schrodinger equation 
is solved numerically ||4] starting with the boundary condition at r = b and ending at a 
value beyond the range of the nuclear potential. At this point the S matrix is obtained 
by matching this numerical solution $ to the required asymptotic form 



$ as (r) ~ 



H 2 (pr) + H^S" 



(6) 



Since the Coulomb potential has infinite range one has to match to Coulomb wave functions, 
so Hi and H 2 are diagonal matrices with entries H^\r] a , p a r) p a r) , the Coulomb 

analogues of the spherical Hankel functions. r] a is the relativistic Coulomb parameter 



Va 



«Mab = a'm a /p a . 



(7) 



Written in terms of the standard regular and irregular Coulomb wave functions, Fg(r],pr) 
and Ge(rj,pr), these Hankel functions are defined as 

Ht (Va,Par) = Ft(r) a ,p a r) -iGt{r]a,p a r) , Hf\r] a ,p a r) = F e (r] a ,p a r) +iGt{r) a ,p a r) . (8) 
These Coulomb wave functions read asymptotically 

F e (r) a ,p a r) 
Ge(r) a ,p a r) ' ~™ cos 
where the Coulomb phase shifts <r^ a are given by 

a La = arg T(£+ 1 + irja) 



r—>oo 

~ sin 



,7T 



„7T 



(9) 
(10) 



(11) 

In case the Coulomb interaction is absent in a channel (like nn), one can simply put r] a = 0. 
Then the Coulomb wave functions become 



Fe{0,p) = pje(p) , Gi(0,p) = -pm(p) , 
HP{^p) = phf\p) , Hf\^p) = phf\p) 



(12) 
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in terms of ordinary spherical Bessel, Neumann, and Hankel functions. The matching pro- 
cedure at Too works as follows. The multichannel Wronskian is defined by 



W($i,$ 5 



m 



Hi 



(13) 



where the prime denotes differentiation with respect to r and the "T" denotes transposition 
in channel space. The Wronskian of two arbitrary solutions $1 and $2 is independent of r 
and equal to because of the boundary condition $(0) = 0. Demanding that 



00 j 1 ^ as 



>oo)) 







we obtain for the partial-wave S matrix 

1 „ ,t 



3" 



IT 



imp 



-_H X - $ J 



-1 r 



IT 



-Ho - $ J 



Imp 



-H' 2 



(14) 



(15) 



where the prime on the Hankel functions denotes differentiation with respect to the argu- 
ment pr. Since the matching is to Coulomb wave functions, what is actually obtained in 
this manner is the partial-wave S matrix with respect to the Coulomb force, denoted by 



Sq +mm+n . How to calculate the scattering amplitude is explained in Section V. The final 
step, the calculation of all the observables from the scattering amplitude, is standard |45|-449 
See in particular Ref. [5C] for the case of antinucleon-nucleon scattering. 




III. THE SHORT-RANGE INTERACTION 

The short-range dynamics is treated with the help of a boundary condition applied at a 
distance r = b. This boundary condition, called the P matrix, is the logarithmic derivative 
of the solution matrix $(r) at a distance r = b from the origin 

(16) 

r=b 

The factor b is included in this definition in order to make the P matrix dimensionless. 
The boundary-condition approach to strong interactions goes back to the work of Feshbach 
and Lomon |yj and earlier. The term P matrix was introduced by Jaffe and Low |52|j53|] 
for use with the bag model where at the energies of the eigenstates of the confined quark 
and gluon degrees of freedom the P matrix exhibits poles that are not necessarily also 
present in the S matrix. (For a review, see Ref. ||54|| .) At this stage, we will not attempt to 
make any connection with multiquark states. In general, this boundary relates the inner- to 
the outer-region physics. In NN and NN scattering the short-range interaction is essentially 
unknown and has to be treated phenomenologically. The long-range physics one understands 
theoretically much better. The P-matrix formalism provides a useful separation between 
these two regions and has become a powerful tool in analyzing scattering data. In the 
Nijmegen partial- wave analyses of NN scattering data [0-0] the P-matrix method has already 
proven its power E5" . 



For the parametrization of the partial-wave P matrix, it is convenient to use a very simple 
model for the short-range dynamics. We assume that the interaction in each partial wave 
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can be described by a spherical well, a potential which may depend on spin, isospin, and 
energy, but which is constant as a function of distance. The P matrix for such a potential 
can be evaluated analytically. For a single- channel spherical-well problem in a partial wave 
with orbital angular momentum i it is given by 

P t = p'bJ' e (p'b)/J e (p'b) , (17) 

where Je{p) = pje(p) with je(p) the spherical Bessel function and p' 2 = p 2 — 2mV, V being 
the depth of the spherical-well potential. The prime on the Bessel function denotes differenti- 
ation with respect to the argument. In the case of NN scattering we take these short-range 
potentials to be complex in order to account for the annihilation into mesonic channels. 
The short-range interaction is in this manner described by a simple state-dependent optical 



potential lp6fl . It turns out that we can take the short-range spherical- well potential inde- 
pendent of the energy and still fit the data properly. This was not possible in the Nijmegen 
partial-wave analyses of the NN data: there we had to make the short-range potential energy 
dependent. That this is not necessary for NN scattering is probably because of the absence 
of high-quality data at low energies such as present in NN scattering. It is, however, crucial 
that we take the real parts to be dependent on spin and isospin. It also turned out that the 
imaginary parts of these short-range potentials could be taken independent of isospin. Each 
partial wave is thus parametrized by a maximum of 3 parameters, a complex spherical well 
for both isospin and 1, where the imaginary parts are equal for both isospins. How many 
and which parameters are actually needed in each individual partial wave is discussed below 
in Section IX. 

In the Nijmegen analyses of NN scattering data, a value of b = 1.4 fm for the boundary 
radius was found to be suitable. In the NN case, the results are rather sensitive on the 
choice of the value of b. The best results are obtained with b = 1.3 fm. Since we use in the 
outer region a real potential consisting of an electromagnetic and a meson-exchange part, 
the coupling to the mesonic channels is completely absorbed in the boundary condition. The 
radius b is therefore a clear measure for the range of the annihilation potential. The fact 
that the results are quite sensitive on b shows that this range is in fact approximately 1.3 
fm. So we find definite indication that annihilation in pp scattering is a rather long-range 
process. 

It remains to discuss the parametrization of the P matrix for the states with £ = J ± 1 
coupled by the tensor force. In these cases it is convenient to use the method also employed 
in the NN case. For a certain value of the isospin, we start with a diagonal 2 x 2 P matrix 
and use additional parameters to describe the short-range mixing between the two coupled 
partial waves, as follows 



cos 9 sin 9 \ I Pi W cos 9 — sin 9 
■ sin 9 cos 9 I \ Po ) \ sin 9 cos 9 



(18) 



The mixing angle 9 can, if necessary, be parametrized as a function of energy, but again in 
our analysis we can take it in all cases independent of the energy. We need these mixing 
angles only for the isospin 1 = states, for the 1=1 states they can be set to zero. 

Time-reversal invariance allows us to choose the phases of the physical states in such 
a way that the potential matrix is symmetric. The full S matrix, including all mesonic 
channels, is then also symmetric and of course unitary. When we restrict ourselves to the 
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NN channels, then this sub-block of the S matrix is of course still symmetric, but no longer 
unitary. This reflects the disappearance of probability (flux) into the mesonic channels. 
Correspondingly, the P matrix in our case is still symmetric, but not hermitian, as in NN 
scattering (below the pion-production threshold). 

Why is the P-matrix approach to partial-wave analyses so convenient and powerful? One 
of the main reasons is that it allows an easy parametrization of the energy dependence of 
the scattering amplitudes. This can be seen as follows. The long-range interactions lead 
to rapid variations with energy of the amplitudes. These variations are much easier to 
parametrize when one uses the P matrix than the S matrix (or K matrix). For instance, 
in the presence of the Coulomb interaction, the S matrix has an essential singularity and a 
branchpoint at zero energy. However, if the Coulomb potential is included in the potential 
tail, these singularities and the corresponding left-hand cut are absent from the P matrix. 
Similarly, the left-hand cuts due to all meson exchanges included in the potential tail are 
absent from the P matrix. All these cuts, however, are present in the S matrix, in addition 
to the kinematical unitarity cut. The point is that the slow variations with energy of the 
amplitudes due to short-range interactions are easy to parametrize, once the rapid variations 
have been taken care of by explicitly including the corresponding long-range potentials in 
the Schrodinger equation. Of course, some left-hand cuts remain, such as the cut due to 
uncorrelated two-pion exchange, as well as right-hand cuts due to the coupling to inelastic 
channels. These, however, lead to much slower energy variations of the amplitudes than the 
long-range electromagnetic interactions and one-pion exchange. 



IV. THE POTENTIAL TAIL 

In order to obtain a good description of the data, the long-range interaction between the 
particles, consisting of the electromagnetic interaction and the one-pion-exchange potential, 
must be included properly. These interactions are model independent in the sense that they 
are (or at least should be) the same in all models of the (anti)nucleon-nucleon force. The fit 
to the data is improved when realistic meson-exchange forces are used as intermediate-range 
potential. 

The long-range electromagnetic potential is included to order a = e 2 /4n, the fine- 
structure constant. The one-photon-exchange potential is derived from the phenomeno- 
logical electromagnetic Lagrangian 

£ 7 = e Q [hf^] A" + e ^ tfa^] - PA") , (19) 

where Q is the nucleon charge in units of e > and k is the anomalous magnetic moment, 
for proton and neutron /i p — 1 + k p — 2.793, and ji n = K n = —1.913, respectively, ip is the 
proton or neutron field and A^ the photon field. If the spatial extension of the nucleons 
is taken into account, the charge and magnetic moments are in momentum space replaced 
by the Dirac form factor Fi(t) and the Pauli form factor F 2 (t), which are functions of the 
four-momentum transfer t. The static limits t — are 

Ff (0) = 1 , i?(0) = , Ff (0) = , F 2 "(0) = . (20) 
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In the point-particle approximation the momentum dependence of these form factors, re- 
flecting the inner structure of the nucleons, is neglected. In this approximation we end up 
with the spin-dependent one-photon-exchange potentials for r > b 



a 



V 7 (r) = -- + 



4M p 2 



a 9 

~Z <J12 



+ 



8/ip — 2 a 



4M p 2 



L • S for pp 



pp 



and 



a 



4Af 2 r 3 



S12 for nn 



nn 



(21) 



(22) 



These potentials are obtained by calculating the one-photon-exchange diagrams in momen- 
tum space and applying a Fourier transformation to configuration space. The momentum 
dependence of the form factors can be taken into account. This leads to short-range mod- 
ifications of the one-photon-exchange potential. Since we use only the potential outside 
r = b = 1.3 fm, we did not include these effects. The use of a' in the central potential for 
pp — > pp takes care of the main relativistic corrections to the Coulomb potential [|57i|58f . It 
is given by 



a 



a 2p/(Mv Uh ) 



(23) 



where t>i a b is the velocity of the antiproton in the laboratory system. In order to appreciate 
the order of magnitude of this factor, at 600 MeV/c for instance, where t>i a b = 0.54, one has 
a' = 1.135 a. The spin-orbit potential comes from the interaction of the magnetic moment 
of one particle with the Coulomb field of the other particle (and is consequently absent in 
nn — > nn). It includes a relativistic correction due to the Thomas precession. The tensor 
potential comes from the interaction of the two magnetic moments. In our energy range the 
Coulomb and the magnetic-moment interaction are the dominant electromagnetic effects. 
The vacuum-polarization potential |59| , which is important |J in low-energy pp scattering 
(below Ti a b = 30 MeV), has a negligible influence here. Two-photon-exchange effects [58 
are neglected as well. 

The following simple one-pion-exchange potential without a form factor is used 



K(r) = /. 



M 



m 



NN 



mr 



+ 



mr 



(24) 



The mass difference between the neutral 7r° and charged ir^ pion is included, so we take 
m = m n o for the elastic reactions pp — > pp and nn — > nn, and m = m n ± for the off-diagonal 
charge-exchange transitions pp <-> nn. In principle, the pion-nucleon coupling constant is 
charge dependent due to the mass difference between the up and down quarks and due to 
the electromagnetic corrections. We introduce here the relevant coupling constants f NNn at 
the different vertices 



/p — /pp7r° ' f n — i nn7r° ' ^fc — /pn7r + /np-;T 



(25) 



for the reactions pp — > pp, nn — > nn, and pp <-> nn respectively. In the PWA a charge- 
independent pion-nucleon coupling constant is used, where fp=f^ = / c 2 = 0.0745, taken 
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from the most recent Nijmegen pp PWA [jBOfl . In Sect. IX, however, this assumption is 
relaxed when we determine the charged-pion coupling constant from the charge-exchange 
data by adding / c 2 free parameter. 

To improve the description of the data, some sort of heavy-meson-exchange potential 
has to be added to one-pion exchange in the outer region. Since the tails of realistic NN 
potentials are remarkably similar, it probably does not matter very much which potential 
one picks, provided it gives a good description of the NN scattering data. We have opted 
for the charge-conjugation-transformed version of the Nijmegen one-boson-exchange soft- 
core NN potential [3D], which is one of the best NN potentials available. (The use of C 
conjugation rather than G conjugation is more natural when one works on the physical 
particle basis.) The tail of this potential has already been used in the Nijmegen pp and np 
partial-wave analyses, and the tail of the corresponding Nijmegen soft-core hyperon-nucleon 
potential [0,0] has been used by us in our PWA of the reaction pp — > AA |3T| , |32|] . The 
following heavy-boson-exchange potentials are included. 

• Pseudoscalar-meson exchange. Included are 7/(549) and r/(958) exchange. For the 
pseudoscalars, including the pion, we use the pseudovector type of Lagrangian 



C 



PV 



47T- 



(26) 



Although equivalent to the pseudoscalar type of interaction for one-meson exchange between 
protons, the pseudovector interaction is favored because it gives a more reasonable two-pion- 
exchange potential and because it leads to at most small breakings of flavor symmetry for 
the coupling constants of the pseudoscalar nonet to baryons |3J] . The scaling mass m n ± is 
conventionally introduced to make the coupling constant / dimensionless. 

• Vector-meson exchange. Included are p(770), u(782), and 0(1019) exchange. The 
Lagrangian (in terms of rationalized coupling constants) is 



£ v = V^g K + ^7^- GK^l (^ $ v - P&v) ■ 



(27) 



• Scalar-meson exchange. We include a (783), / (975), and /q(760) exchange. The 
Lagrangian is 



(28) 



The scalar mesons have always been a controversial topic. In early one-boson-exchange 
models for the NN interaction there was a clear need for an isoscalar scalar "<r" meson 
with an effective mass of about 550 MeV ]|63H 65|]. While no such low-mass particle exists, 



there was some evidence in production experiments for a broad structure e(760) under 
the p°, often explained away as a strong tttt final-state interaction. Later it was pointed 



out 



that such a wide (r ~ 640 MeV) e(7Q0) simulates the narrow- "cr" exchange 
in one-boson-exchange models for the AW interaction. 

The situation in phase-shift analyses of tttt scattering data, obtained from reactions as 
ttN — > 7T7rA, has for a long time been confusing and not conclusive [p8|] . In these analyses 



the assumption has always been that only tc exchange is relevant, while a\ exchange can 
be neglected. Very recently, however, the situation has been much clarified ||69|| . Data on 
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ttN f — > 7t + tt~N using a polarized target provide unambiguous evidence for a broad 7 = 
++ (750) state, when a proper amplitude analysis is done, including also exchange. In a 
similar amplitude analysis of data on K + n f— > K +r K~p [70| evidence is found for 1 = 1/2 
+ (887) strange scalar mesons under the if* (892). 

In the quark model, several mechanisms give rise to scalar (J p = + ) mesons. The 
simplest model is the 3 Po QQ states. Then there are the glueball states and the cryptoexotic 

2 2 

Q Q 2 states ||71|| . A physical scalar meson will in general be a mixture of QQ, Q Q 2 , and 
glueball components. The QQ states are expected JF^] near the other 3 P QQ mesons, that 
is around 1250 MeV. Glueballs are also not very likely to exist below 1000 MeV |73||. For 

t2. 



the Q Q 2 states, however, one [[H]] does predict a low-lying nonet of scalar mesons. The 
lowest state, with only nonstrange quarks, has 1 = and decays into tttt. It can be identified 
with the £(760) under the p°. This nonet contains also a nearly degenerate set of I = and 
1=1 cryptoexotic scalar mesons (like the p(770) and a; (782)) with an ss pair. These are 
easily identified as the / (975) and a (983) mesons, previously called S'*(975) and 5(983) 
respectively, with their relatively large branching ratios into if if . The nonet is completed 
by a set of broad I = 1/2 strange mesons if,* (887) seen [74] under the if* (892). 

• Next to these conventional mesons, the Nijmegen soft-core potential (originally derived 
from Regge-pole theory) also contains pomeron exchange, which in QCD is understood as 
color-singlet two- or multigluon exchange J75|-77], and the weak diffractive scalar-like part of 
the tensor-meson exchanges. The short-range non-local terms of the potential are neglected, 
which is a very good approximation outside r = b = 1.3 fm. 

The consequences of meson exchanges for the AW interaction have been examined by 



many authors, for instance by Dover and Richard ITq-pOl. It turns out from a qualitative 



investigation that while in the AW case there is a strong coherence between the isospin 
1=1 spin-orbit forces, in the NN case very strong tensor forces occur, especially for isospin 
1 = 0. In the NN system the vector a;(782) and the scalar /q(760) exchange make up the 
strong spin-orbit force that splits the 3 -Po,i,2 phase shifts, but the central potentials of these 
exchanges largely cancel each other. Similarly, the spin-orbit forces due to the exchange of 
the vector p(770) and the scalar a (983) add up, but the central potentials cancel. Applying 
charge conjugation to the different meson exchanges, one sees that in the NN potential 
the central potentials from u;(782) and /q(760) exchange add coherently to a very strong 
attractive potential. This has led to speculation about the existence of AW bound states 
and resonances |8~T|-|S5|j . The tensor potentials from p(770) and 7r(140) exchange also add up 
and dominate the charge-exchange reaction pp — > nn. The same phenomenon is present in 
the reaction pp — > AA, where if (494) and if* (892) exchange conspire to build up the strong 
tensor force that is the hallmark of this reaction ||31|| . 



V. THE AMPLITUDES 

The evaluation of the scattering amplitude in the presence of electromagnetic effects is 
a non-trivial problem that requires special care. In the pp PWA the precise values and 
energy dependence of the phase shifts are influenced significantly by the electromagnetic 
interaction. The same will obviously also be true in the pp case. 

We start by discussing the case where next to the nuclear force we only have the Coulomb 
potential to deal with. The scattering amplitude due to the Coulomb force is given by 
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(s'm'a , \M c (9)\sma) = -^^^ - Va e -^h^e )+ ^ , a 

p a {l ~ COs9) 

= -5 ss '5 mm ,5 aal ^- a{sin 2i_ d)1+lVa ■ ( 29 ) 

The partial- wave decomposition of Mc(9) in terms of Coulomb phase shifts and Legendre 
polynomials does not converge point-like, due to the infinite range of the Coulomb potential. 
However, it can be summed in the sense of distributions |p6|-p8|] to give the Coulomb ampli- 
tude Eqn. (P5p. In order to make a partial- wave decomposition of the scattering amplitude, 
the total scattering amplitude M c +n(9) is split as follows 

M c+N {9) = M c (6) + Mg +N (6) , (30) 

where the amplitude Mq +n {6) is the nuclear scattering amplitude in the presence of the 
Coulomb potential. Its partial-wave decomposition reads 



(s'm'a'\MS +N (8)\sma) = £ ^4n(2£ + 1) & e C^J m Cl_ m , * J m Y*_ m ,(9) 

ti'j 

(I's'a'lS^iSS+N - l)S 1 J 2 \£sa)/2tp a , (31) 

where Sq +n is the nuclear S matrix in the presence of the Coulomb force. The total S 
matrix Sq+n, due to the Coulomb and nuclear interaction, is then given by 

c _ c 1 / 2 qC a 1 / 2 (on\ 
•JC+N — iJc ^C'+N D C* • 

Sc is the Coulomb S matrix with matrix elements 

(£'s'a'\S c \£sa) = 5 u >5 ss >5 aa > exp(2ia^ a ) . (33) 

The nuclear S matrix in the presence of the Coulomb force Sq +n is obtained by solving the 
Schrodinger equation numerically and matching to Coulomb wave functions (see Section II). 
It is given by Eqn. (|15|). 

Next we discuss the generalization to the case where next to the Coulomb force also the 
magnetic-moment interaction is present . Although this latter potential has a finite range 



and consequently the partial-wave decomposition converges, it is still much more practical 
to split off in the scattering amplitude the contribution of the magnetic-moment interaction. 
In this way the magnetic-moment interaction can be included in all partial waves and the 
summation in the nuclear amplitude converges much faster. We thus write 

M C+MM+N {9) = M c (6) + MS +MM (9) + MSt^ +N (9) . (34) 

Here Mq +mm (9) is the scattering amplitude of the magnetic- moment interaction in the 
presence of the Coulomb force. Since this amplitude is almost exactly in phase with the 
Coulomb scattering amplitude Mc{9), it is essential that the effect of the magnetic-moment 
interaction is evaluated in Coulomb-distorted- wave Born approximation (CDWBA), and not 
in plane- wave Born approximation, as was pointed out by Knutson and Chiang |p0[ . The 
contribution, in CDWBA, of the magnetic-moment interaction to the K matrix is 
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M a f°° 

s'a'\K MM \£sa) = -5 aa > — / dr F ef (ri a/ ,p a >r)V M M{r)Fi(r) a ,p a r) . (35) 

Pa Jo 



Integrals of this type can be evaluated very rapidly and accurately by a backward-recursion 
algorithm ||91{| . The magnetic-moment interaction Vmm is given in Eqn. fl2"T| ) for pp —>■ pp and 
in Eqn. fl22|) for nn — > rm. All these partial-wave contributions are subsequently summed to 
obtain the total amplitude M^ +MM (9). In practice it turns out that the spin-orbit poten- 
tial of the magnetic-moment interaction leads to a contribution Z^s to (lla\M^ +N (6)\10a) 
that converges much too slowly to be summed term by term. This part can be summed 
analytically ptJ. For antiproton-proton scattering the result is 



J LS 



where we have defined 



M P hs f ^la-cose) _ i (1 _ cos ^ ; (36) 



sin^v^V 2 



hs = -T^( 8 ^- 2 ) • (37) 
There is a similar contribution —Zls to (10a\M^j +N (6)\lla) . The corresponding result for 



pp scattering |90(] is the properly symmetrized version of this expression. The partial-wave 
decomposition of M c +mm+n(@)i the nuclear scattering amplitude in the presence of the 
Coulomb force and the magnetic- moment interaction, is similar to Eqn. (pi]), but there now 
appears the following S matrix 

C _ C-Va/oC A 1 / 2 qC+MM (qC ql/2 /nM 

&C+MM+N — >-> c yo c+MM j o c+MM+N yoc+MM J °c ■ \°°) 

Since the magnetic- moment interaction contains a tensor part, the matrix Sq +mm is not 
diagonal in orbital angular momentum. However, the square root of this matrix is still 
well-defined. What it all comes down to, is to rewrite the S matrix in the following manner 
in order to split off the Coulomb amplitude Mc and the amplitude Mc+mm due to the 
magnetic-moment interaction in the presence of the Coulomb force 

Sc+MM+N — 1 = (Sc — 1) + S l J 2 (Sc+MM ~ l) S(J 2 + 

C 1 / 2 ( qC \ l / 2 (qC+MM AfqC V /2 qV2 /oq\ 

°C \ D C+MMJ \°C+MM+N l )\ D C+MMJ °C j 

Other electromagnetic effects, like vacuum polarization, can be treated in a similar way. For 
a more extensive discussion we refer to Ref. ||89|| . Because of the long range of the magnetic- 



moment interaction the matrix elements of Sq+mm+n are hard to calculate. We will use the 
approximation 

cC+MM qC r m \ 

°C+MM+N ~ °C+N ■ l^ u J 

Exactly the same approximation is made in the Nijmegen AW partial-wave analyses |89 . 
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VI. THE PHASE SHIFTS 



In this section we give the parametrization of the NN S matrix in terms of phase-shift 
and inelasticity parameters. We first seek guidance in the way this is done in analyses of NN 
scattering below the pion-production threshold where the S matrix is unitary and symmetric. 
The symmetry of the S matrix is a consequence of time-reversal invariance which allows one 
to choose the phases of the in- and out-states such that the coupled-channels potential 
matrix, and thus the S matrix, is symmetric. If there is conservation of flux, the S matrix 
is unitary. 

The phase shift for uncoupled partial waves with £ = J, s = 0, 1 is defined by parametriz- 
ing the 1 x 1 S matrix as 

S J = exp(2i5) . (41) 

One usually denotes the different phase shifts by 5i for singlet s = waves, and by 5e,j for 
triplet s = 1 waves. For the partial waves with £ = J±l, s = 1 coupled by a tensor force one 
writes the 2 x 2 S matrix in terms of two phase shifts £ j+i,j and one mixing parameter 



Ej. A popular parametrization is the "eigenphase" convention of Blatt and Biedenharn |92 
in which the symmetric S matrix is diagonalized by way of a rotation 

S J = exp(-iEjfjy) exp(2i<5) exp(ieja y ) , (42) 

where 

(43) 




More often one uses the "bar-phase" convention of Stapp, Ypsilantis, and Metropolis |[93|| , 
in which 

S J = exp(iS) exp(2iEja x ) exp(i5) . (44) 

An advantage of the "bar-phase" convention is that the parameters go to zero when the 
interaction vanishes, unlike the mixing parameter in the "eigenphase" convention. Only in 
the former case is the mixing parameter a measure of the strength of the off-diagonal tensor 
force. We will use the "bar-phase" parametrization for the three elastic parameters, since 
this is at present the common choice in analyses of NN scattering. 

In the presence of coupling to annihilation channels the S matrix describing NN scat- 
tering is only a submatrix of the much larger coupled-channels 5* matrix, and is therefore 
still symmetric, but no longer unitary. This doubles the number of parameters needed. The 
parametrization of the S matrix for uncoupled partial waves with £ = J, s = 0, 1 in the 
presence of annihilation requires two parameters. One writes 

S J = r] exp(2i5) , (45) 

with |^| < 1. To denote the inelasticities for the different partial waves we will use a similar 
notation as for the phase shifts: rji for singlet s = waves, and tj^j for triplet s = 1 waves. 
For the coupled partial waves with I — J±l, s — 1 six parameters are needed, and it is not 
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so easy to think of a convenient parametrization which satisfies all constraints from unitarity, 
is completely general, and free from non-trivial ambiguities. Fortunately, the essential work 
has already been done by Bryan and others [p5"-9~T|| in the case of NN scattering above 



the pion-production channel. Bryan generalizes the "bar-phase" convention by writing 

S J = exp(i5) exp{ieja x ) H J exp{ieja x ) exp(i5) , (46) 

where if is a three-parameter real and symmetric matrix representing inelasticity. Bryan 
calls this matrix N, but we use H (capital rj) to stress the analogy with the case of uncoupled 
partial waves written in the form 

S J = exp(iS) 7] exp(iS) . (47) 

If the inelasticity vanishes, then H tends to the unit matrix, and the "bar-phase" 
parametrization is recovered. There are several nice ways to parametrize H, but we find it 
convenient to follow Klarsfeld |56|, who diagonalizes H in Blatt-Biedenharn fashion 



H J = exp(-iujja y ) ( Vj Q 1J ^ ^ ^ exp(iu jC r y ) , (48) 

where the diagonal matrix contains the "eigeninelasticities" f]j-i,j and r]j +1 j. In this way 
the partial-wave annihilation cross section depends only on the "eigeninelasticities," since it 
is proportional to 

Tr (1 - H 2 ) = 2 - V 2 j_ hJ - V 2 J+ltJ . (49) 

If the phase parameters are actually searched for on a computer, it is better to write all 
inelasticities as 77 = cos 2p. In this way, all parameters are real and unbounded. The mixing 
parameter ujj is finite as the inelasticity vanishes, just as the mixing parameter ej in the 
"eigenphase" convention tends to a finite value when the off-diagonal tensor force goes to 
zero. Sprung p7 has extended the Bryan parametrization to allow the use of "eigenphases." 



In this case one writes 

S J = exp(—iejcr y ) exp(i5) H J exp(i5) exp(ieja y ) . (50) 



The matrix H in that case is, in general, not equal to the matrix H in Eqn. Q4"B]). As stated 
above, we will use the "bar-phases." The algorithm, due to Bryan |9~4}1 , to extract the phase 
parameters and inelasticities from the S matrix presented in numerical form is reviewed in 
the Appendix. 



VII. STATISTICS 

Statistics is an essential ingredient in analyses of large amounts of scattering data. The 
theoretical predictions are compared with the experimental data using a least-squares fitting 
procedure in which the model parameters are adjusted to the data. During this process one 
continually scrutinizes the data and passes sentence on the quality of different sets, which 
sometimes have to be rejected on the basis of statistical criteria. This goes on until a 
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final verdict is reached and the building of the data set is completed. In this section this 
procedure is outlined and the statistical tools required for our purpose are presented. A 
more exhaustive treatment can be found in Ref. ||. 

We start by assuming for the moment that the measurements have no normalization 
uncertainties and that no other type of systematic error is present. In a certain experiment, 
denoted by a subscript a, one has measured N a data points. Each such measurement with 
statistical error is denoted by E a ^ ± e a ^ {i — 1, . . . , N a ). The model prediction for a certain 
data point is given as M a i (p), where the model parameters are arranged in a vector p, with 
entries p a (a — 1, . . . , N pa _ r ). The parameters are adjusted to the data by minimizing the 
X 2 -function 



x 2 (p) 



Ex'(p) = EE 



a i=l 



M a>i (p) - E a , 



(51) 



with respect to all parameters. 

In practice, however, measurements usually do have an overall normalization uncertainty, 
specified by the experimentalists. These errors can be taken care of by introducing for each 
group of data a normalization parameter v a with error e a $. The normalization of a certain 
group is then given as v a = 1 ± e a> o . This means essentially that for each group with a finite 
normalization uncertainty another free parameter v a has been introduced, to be determined 
in the fit, as well as an additional datum 1 with error e a0 . Since we want to restrict the 
parameter space to the model parameters only, we employ the following trick to take the 
normalization parameters into account. We redefine the x 2 -function as follows 



X 



(P) = EXa(P) 



N a 

E min E 



- 2 




+ 


£a,0 



1 2 



(52) 



In this way, the normalizations are adjusted trivially, by minimizing in each iteration a 
quadratic function. In case the normalization of an experiment is completely unknown, 
e a fi = oo and we can remove the second term on the right-hand side of Eqn. (|52"D. The cor- 
responding normalization v a is determined in the fit. We say the normalization is "floated." 
If the normalization is exactly known we fix the normalization at v a = 1 and we can remove 
again the second term. Angle-dependent normalization errors can be treated in a similar 
manner. The minimum value Xmm = X 2 (p) lp=p n ;« ^ s reached when 



d X 2 (p)/d Pa = , (53) 
for all values of a. At this point one defines the error matrix E of the parameters as 



(E 



- X U = 2 9 2 X 2 (p)/d Pa d Pp 



P = Pn 



(54) 



The standard error on the parameter p a is (E aoi ) l l 2 . Assuming that the x 2 "f unc tion is 
quadratic near its minimum, one can show that this is the variation in p a that gives a rise 
Ax 2 = 1 in Xmiri; when the remaining parameters are refitted [98]. 

According to this discussion the following integers can be defined. In the fit one must 
determine A^ n normalization parameters as well as A^ par model parameters. The actual num- 
ber of free parameters is therefore A^ p = A^ par + A^ n . Of these A^ n normalization parameters 
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N ne have a finite normalization error and the rest N n f = N n — N ne is the number of "floated" 
normalizations. The total number of experimental scattering observables is N Q ^ S and the 
actual number of scattering data is iV data = N ohs + N ne . The number of degrees of freedom 
is iV df = iV data - iVfp = W obs - AW - N Qi . 

In the process of screening the database, we employ certain rejection criteria to remove 
data points that spoil the statistical quality of the data set, for instance by underestimated 
statistical errors or by unspecified systematic errors. When statistical errors are underesti- 
mated, data pretend to give more information than they actually do, causing false results. 
By systematic errors we mean those errors that do not average to zero when the measure- 
ment is repeated many times, causing them to have a correlated effect on the data. When 
not treated correctly, systematic errors can bias the values of the model parameters and 
therefore also of the phase parameter. A detailed discussion of this point can again be found 
in Ref. ||. Of course, we only reject data if there is conclusive evidence against them. The 
rejection criteria are constructed in such a way that in a purely statistical ensemble they 
have a very small chance to occur. The following rejection criteria have been employed by 
us. 

(i) An individual data point E a j that has x\% > 9 is rejected. Rejection of such outliers 
will give more accurate values for the model parameters. This rejection criterium corresponds 
to the usual three-standard-deviation rule. It implies that a correct datum has at most a 
chance of 0.27% to be rejected. The same criterium applies to an experimental normalization. 
If it contributes more than 9 to Xmm> tins datum is rejected and the normalization is floated. 

(ii) A group of iV a data is rejected if its x\ i n the multienergy fit is less than a minimum 
X\aw(Na)- The values for X\ ow can be found in Ref. ||. They are constructed such that again 
the chance for a correct group to be rejected is at most 0.27%. This rejection criterium is 
used as a means to avoid systematic errors. Very low values of x 2 are very probably caused 
by systematic errors present in the data and are not a virtue of a model. This criterium 
generally does not seem to be appreciated by the uninitiated in statistical methods. 

(Hi) A group of A" a data is also rejected if its x\ lIy ^ ne multienergy fit exceeds a maximum 
Pxl lgh (N a ), where (3 = N di /(N di + N pav ). The values for 

Xhigh can a l so be found in Ref. ||. 

VIII. THE NIJMEGEN 1993 ANTIPROTON-PROTON DATABASE 

A. Set-up of the 1993 database 

In order to perform a pp PWA, we first had to put together a statistically sound data 
set. In NN analyses, one has over 30 years experience with the data, and by now a proper 
database is more-or-less agreed upon. In the pp case, we must start from scratch, but 
fortunately we can use the experience of the NN analyses to arrive at an acceptable pp 
database. Let us summarize the general features of our database. We will include in our pp 
database all available pp scattering data below antiproton laboratory momentum pi a b = 925 
MeV/c published in a regular physics journal since 1968. We do not take into account data 
published only in conference abstracts, in conference proceedings, and/or theses. The reason 
we restrict ourselves to this momentum range is that it corresponds roughly to the energy 
range of the AW partial-wave analyses below kinetic energy Ti a b = 350 MeV (momentum 
925 MeV/c corresponds to Ti a b = 379 MeV) and that we want to include the accurate 
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backward elastic cross sections between 406 and 922 MeV/c taken by Alston-Garnjost et 
al. ||117|| . Low-energy data, below p lab = 175 MeV/c, so T lab ps 15 MeV, are of course 
almost nonexistent here. Only scattering observables are analyzed, other "data" like for 
instance the real-to-imaginary ratio of the forward scattering amplitude or the slope of the 
pp forward nuclear amplitude are omitted, since the extraction of these quantities from the 
data is model dependent. A summary of the Nijmegen 1993 pp database can be found in 
Table § 

Several experiments ||1 33| |1 35f| have reported a resonant structure in the antiproton- 
proton total cross section near 490 MeV/c without agreeing, however, on the exact position, 
strength, and width of a possible resonance. We do not include these data in the PWA, 
since more recent and accurate measurements of total cross sections ||1 3 6| . |1 1 1| , |1 9|J1 3 7| , |1 1 4 



have convincingly ruled out the existence of this type of rather broad resonance. (A more 
narrow type of resonance, however, is not ruled out by the existing data. In fact, there is 
some statistical evidence for a narrow structure in backward elastic cross sections around 
509 MeV/c [|138|| . Probably, only an accurate measurement of the backward charge-exchange 
cross section can definitely settle this issue.) 

High-quality total cross sections below 400 MeV/c have been measured at LEAR by 
the PS172 collaboration The effects of the pure Coulomb force and Coulomb-nuclear 

interference are not significant, except perhaps at low momenta < 300 MeV/c. This may 
be the reason why we had to reject the data at the two lowest momenta. 

In a number of experiments [|134| , |139| , |111|J140| - |144|| the pp annihilation cross section into 
charged mesons has been measured. We cannot include these data, not even with a floated 
normalization, since the momentum dependence of the cross section for annihilation into 
neutral mesons is not known. At LEAR the total annihilation cross section, including 
annihilation into neutral channels, has been measured from 180 to 600 MeV/c by the PS173 
These last data are included. 



105.100 



group 

We do not include integrated elastic-cross-section data | |1 10| . |104| , |145| . |1 13j . |l 15| ] , due to dif- 
ficulties in a proper treatment of Coulomb and Coulomb-nuclear interference. Integrated 
charge-exchange cross sections, on the other hand, are included in the database. In two 
experiments ||108| , |9~9[1 this observable was measured over a large momentum region. For the 
experiment of Ref. |)9| data points at the lowest momenta are rejected since they are in 
conflict with more recent measurements done at LEAR by the PS 173 group ||103|| . Our 
PWA clearly favors this last experiment. 

In the pre-LEAR era a very large part of the experimental data on pp scattering consisted 
of elastic differential cross sections ||1 1CI|J104|J129| , |132| , |T28| , |1 13|| . The most accurate data were 
those taken by Eisenhandler et al. at 690, 790, and 860 MeV/c (and higher momenta) ||128|| . 
With the exception of the data at 194.8 MeV/c of Ref. [|10§] and the data at 910 MeV/c of 
Ref. [|132| we find these pre-LEAR data to be consistent. An accurate measurement of the 
backward elastic differential cross section at cos 6* = —0.994 was done by Alston-Garnjost 
et al. ||117|| between 406 and 922 MeV/c. Since 1983, differential cross sections on pp — > pp 
have been measured by different groups at LEAR and at KEK. For an extensive discussion 
of these data we refer to the next subsection. 

A number of pre-LEAR experiments determined the elastic differential cross sec- 
tion at forward scattering angles in order to study the Coulomb-nuclear-interference re- 
gion p6l[Uq , |iT^ , |l47|1 . At LEAR this was done by the collaborations PS172 PT7|JI25B and 
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PS173 ||148|| . In general, all these data are well described in the PWA, although in some cases 
data points in the extreme-forward region had to be rejected since they are contaminated 
by multiple Coulomb scattering in the target (so-called Moliere scattering). 

Some of the most important results coming from LEAR are the high-quality elas- 
tic analyzing-power A y (polarization) data, measured by the PS172 ||122| , |123|| and 
PS198 [|130| , |12CI|| collaborations at a number of energies above 439 MeV/c. Since in the 
pre-LEAR era only very few and inaccurate data points existed | |132|J127| , |131| , these exper- 
iments mean real progress. We find the data of PS172 and PS198 to be consistent, except 
for their normalizations. In our PWA we float the normalizations of the A y data at 497, 
523, and 679 MeV/c. PS172 has also obtained the first (not very accurate) results on the 
elastic D yy depolarization ||150 |. For these data we do not take a normalization error into 
account, in view of the large error bars. 

Before the advent of LEAR in 1983, also charge-exchange differential observables were 
scarce. Some differential cross sections existed ||1 1 9| , |1 29| , |1 5 1| , |1 2 1| , |1 52j ] , but these were not 
very accurate. Since 1984, however, the situation has improved enormously. High-quality 
differential cross sections have been obtained at KEK between 392 and 781 MeV/c ||116|| , 
and at LEAR by the PS173 group at low momenta between 183 and 590 MeV/c ||103|| . One 
of the most important experiments at LEAR has been PS199 whose goal was to study the 
spin structure of the charge-exchange reaction. So far, it has obtained very accurate data on 
the differential charge-exchange cross section at 693 MeV/c and very important data on the 
charge-exchange analyzing power A y between 546 and 875 MeV/c ||126| , |124|| . More results 



from PS 199 can be expected in the near future. Very recently, an new LEAR experiment 
called PS206 ||153|| has been approved that will further study the charge-exchange reaction. 



B. Flaws in elastic differential cross-section data 



Since 1983 several experiments have measured the elastic differential cross section at 
different momenta. PS173 at LEAR measured this observable at low momenta between 
181 and 590 MeV/c ||101|J102|| . It was subsequently measured at KEK between 392 and 781 
MeV/c P5] , by PS172 at 679, 783, and 886 MeV/c |23j (and at higher energies), and by 
PS198 at 439, 544, and 697 MeV/c pO| , [T20l . Unfortunately, these different experiments 



do not appear to be consistent with each other, nor with the accurate pre-LEAR data of 
Eisenhandler et al. [ 128 | which are described very well in the PWA. In fact, this is the most 
serious flaw in the pp database and a major obstacle in fitting potential models [55]. It is 
highly probable that some of these data contain unspecified systematic errors, or that the 
statistical errors have been underestimated. 

Since we are talking about a total of 540 cross sections (not counting the Eisenhandler 
data) we decided to put some effort in trying to determine what is wrong with these data, 
instead of rejecting them outright. This turned out to be very difficult. The discussion 
about statistical and systematic errors in most of the original papers can be called marginal 
at best. Sometimes results from Legendre fits are presented, but in several cases it turns out 
that the Xmin °f these fits is much too high for a statistical data set. The procedure that we 
followed to examine these data was as follows. We fitted each differential cross section with 
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da/dQ = aiPi(cos9) . 

1=0 

f max was determined by the requirement that the error on the corresponding coefficient 
ci£ was smaller than the coefficient itself. In the ideal situation this would give a fit with 
Xmin/^df ~ 1.0. However, for few groups this was actually the case. We then enlarged the 
errors by adding a point-to-point systematic error in quadrature to the quoted statistical 
errors. The amount of error added was determined by the requirement that now indeed 
Xmin ~ -Wdf- Once for a group this result was approached, outliers were removed (three- 
standard-deviation rule). For some groups special measures had to be taken, which we will 
discuss now. The results of these investigations are summarized in Table |TT]. 



The data measured by PS 173 ||101| , |102|| are at the most forward angles contaminated by 
effects due to multiple scattering in the target (Moliere scattering). Since the experimental- 
ists do not present their data corrected for these effects, and since these corrections depend 
on details of the target used, there is no way for us to take these effects into account. Conse- 
quently, the only sensible thing to do is to reject the data at angles where multiple-scattering 
effects are believed to be seen. These data were used to extract the real-to-imaginary ratio 
of the forward scattering amplitude \\L4S ]. It is perhaps remarkable that in Ref. ||148|1 values 



for this ratio were presented at more momenta than the four for which the corresponding 
differential cross sections were presented in Refs. [|101| , |102j 1. These remaining cross sections 



have never been published. After removing the points at forword angles, difficulties remain. 
At 287 MeV/c it is necessary to reject five individual data points and add a 5% point-to- 
point error. At 505 MeV/c two outliers must be removed and at 590 MeV/c a 3% error 
must be added. 

The data taken at KEK [|115|| present serious difficulties of a different kind. The only 
manner to achieve satisfactory results in a Legendre fit seemed to be to either reject a 
large number of data points or to enlarge all errors by a significant amount (~ 6 — 7%). 
However, improvement could be obtained in the following manner. It turned out that most 



difficulties resided in the data taken in the "one-prong" region ||115|| , i.e. at forward and 
backward angles. The data at intermediate angles, the "two-prong" region, appear to be 
less troublesome. In fact, in Ref. ||115|| different systematic errors are quoted for these two 



regions, although the experimentalists perform Legendre fits to the data as one group. When 
we split the groups at 490, 591, 689, and 781 MeV/c into two different parts, it suffices to 
enlarge the statistical errors by a smaller amount, and in some "two-prong" regions the data 
need no corrections at all. 



The description of the data taken by PS172 [ |123j| improves much when the points at the 
most forward and at the most backward angle at removed at all momenta (F. Bradamante, 
private communication). At 886 MeV/c three additional outliers have to be rejected. Espe- 
cially for these groups there remains a big problem with the normalization ||123|| (see below). 
For the data at 697 MeV/c [|130|| from the PS198 group the statistical errors may have been 
underestimated. We had to enlarge these by adding a 6% point-to-point error. 

After all these corrections have been applied to the different groups the next step is to 
see to what extent the Legendre coefficients at comparable momenta are consistent. These 
coefficients are presented in Table |TJ. It is immediately clear that significant problems occur 
in the normalization a® of the different groups (which is why we present the "renormalized" 
Legendre coefficients a^/ao for I > 0). This is especially true for the PS172 data. The 
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probable reason [|123|| for the difficulties is that data could only be taken in a very limited 
angular region, so that only a small fraction « 5% of the cross section is detected. Properly 
normalizing the data is then especially difficult. These data would have to be treated 
with a "floated" normalization. Certainly a 10% normalization uncertainty as suggested in 
Ref. ip.23 ] is not sufficient. Since also the KEK data in the "two-prong" region cover a limited 
angular region it would be a good idea to float these normalizations as well, since they are 
not consistent with the normalizations of the data in the "one-prong" regions. Apart from 
these normalization problems it is clear from a study of Table |Tl| that it is very unlikely 
that the different experiments are consistent. Certainly, the PS173 data at 590 MeV/c are 
not compatible with the KEK data at 591 MeV/c; the PS172 data at 783 MeV/c are not 
compatible with the KEK data at 781 MeV/c, nor with the Eisenhandler data at 790 MeV/c. 
The PS198 data at 439 MeV/c do not appear to intrapolate between the KEK data at 392 
and 490 MeV/c, etcetera. So, although many beautiful data have come out of LEAR and 
KEK, unfortunately elastic differential cross sections are not among them. 

To summarize, at present we are prejudiced in favor of the pre-LEAR data by Eisen- 
handler et al. ||128|| , although we cannot exclude completely the possibility that these data 
contain unspecified systematic errors and that one of the LEAR or (more unlikely) the KEK 
experiments is correct after all. Preliminary study showed that we could obtain a reasonable 
fit to the LEAR data from PS172 and PS198, but only at the cost of rejecting the pre-LEAR 



data of Eisenhandler et al. ||128|| and of Sakamoto et al. ||113 |. And still, the problems with 



the data from PS173 and KEK would remain. Further investigation is required before such 
a drastic step will be taken. Obviously, a dedicated new experiment that might shed some 
light on this issue would be highly welcome. 



IX. RESULTS 

After deciding on the final content of the Nijmegen 1993 pp database to be used in the 
PWA, the free P-matrix parameters are fitted to these dataQ In Table | we present an 
overview of the results of this fit for all the pp scattering data. In this Table one can find 
the values for Xmm f° r eacn individual group, the normalization predicted by the PWA, and 
the data points and groups that have been rejected by us on statistical grounds. The total 
number of scattering observables that are rejected is 204, not including the problematic 
LEAR and KEK elastic differential cross sections discussed in the previous Section. Three 
normalization data are rejected. We reject three groups because of an improbable high Xmin' 
and two groups because of an improbable low Xmin- m the fi na l fit we have (in the notation 
of Sect. VII) 

Aobs = 3543 , N n = 113 , N ne = 103 , iV par = 30 , 

so that 



1 For the sake of completeness we mention that at the time of our final fit, the data of 
Refs. [151,152,147] were not available to us. Also, we were at that time not aware of the exis- 
tence of the data of Refs. |136.137|. 
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iV data = 3646, W fp = 143, N n{ = 10 , N# = 3503 



When the data set is a perfect statistical ensemble and when the model is totally correct 
one expects 

(xiJN d{ ) = 1.000 ± 0.024 . 
Our best fit to the final 1993 database results in 

xLn = 3801.0 , 

corresponding to 

X mi n/^data = 1-043 and xLJNx = 1.085 . 

The 3543 scattering observables contribute 3700.9 to x min , which means that the 103 nor- 
malization data contribute the remaining 100.1 of Xmin- 

In order to get a feeling for some of these numbers, let us compare them to the AW 
case. In the Nijmegen AW PWA the database contains A*d ata = 4301 AW scattering data, 
which are described with xLm = 4263.8 or x m in/-^data = 0.991 0. The number of model 
parameters needed in our case is 30, which is a reasonable number, in view of the fact that 
21 parameters were used in the Nijmegen pp PWA and an additional 18 in the np PWA. 

The values for the P-matrix parameters and their errors are tabulated in Table fV[ The 
parameters for the higher partial waves are the same as the corresponding state given in this 
Table. For instance, the 1 F% wave has the same parameters as the 1 D2 wave, the wave 
the same as the 3 D2 wave, the the same as the 3 P2, and so on. This is just a convenient 
prescription. For higher partial waves the short-range parametrization is irrelevant due to 
the centrifugal barrier and the dynamics is completely determined by the potential tail. As 
explained in Section III, these parameters correspond to a spin-dependent optical potential 
for the short-range interaction. It can be seen that even if one introduces for each partial 
wave of specific isospin a simple complex spherical-well potential, then still by no means all 
of these parameters can be determined from the existing data. This is in striking contrast 
with for instance the pp analysis where one needs 4 parameters already to describe the 1 Sq 
channel. The reason is that the pp 1 Sq phase shift is very accurately known at very low 
energies (below 3 MeV). A proper description of the P waves in the pp case also requires 
more than one parameter for each wave. In our pp case, no single partial wave (of specific 
isospin) needs more than one parameter for the real part of the short-range interaction, 
but on the other hand much more partial waves contribute significantly to the scattering 
process. This is especially true for the charge-exchange reaction, as can be seen from Table 
where we give the partial-wave cross sections for the elastic and charge-exchange reaction 
at a number of momenta, as well as the total and annihilation cross section. Compared 
to nucleon-nucleon scattering, one notes a large contribution of P and D waves to the 
cross sections already at low momenta. The reason for this is the greater strength of the 
antinucleon-nucleon potentials, especially the central and tensor potentials. 

In Figure [I] we show the results for the total and annihilation cross sections as a function 
of momentum. The high-quality data are from the LEAR experiments PS172 ||106| , |114|| (total 



cross sections) and PS173 [p.05| , |100|1 (annihilation cross sections). We calculate total cross 
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sections with the optical theorem. Obviously, these cross sections can only be compared 
to the experimental data, when the effects of the Coulomb interaction can be neglected. 
Except for the lowest momenta (pi a b < 200 MeV/c) this is probably a good approximation. 
An example of the fit to the differential cross sections from Eisenhandler et al. ||128|| is shown 
in Figure |2|. In Figure [| the fit to the backward elastic cross sections (cos8 = —0.994) from 



Ref. [117 1 is demonstrated. Assuming that these data are correct, it can be seen from this 



Figure that there appears to be room for improvement at the lowest momenta. This is 
precisely the momentum region (pi a b ~ 509 MeV/c) where some statistical evidence for a 
resonance was found []138|| . The differential cross section for charge-exchange scattering is 
given in Figure |], compared to recent high-quality data taken by PS199 [ |126|| at 693 MeV/c. 
This is one of the most constaining experiments in the database. In order to fit this group 
properly, orbital angular momenta up to i = 10 must be taken into account. The cross 
section exhibits the typical dip-bump structure at forward angles, which can be understood 
as an interference effect between one-pion exchange and a background ||1 54|| . Unfortunately, 
no data of similar quality have been taken in this dip-bump region. 

In Figure [5] we give the results for the analyzing power (polarization) for elastic scattering, 
compared with the recent data from PS172 |22| and from PS198 [p0l , |T20] . The fits to the 
analyzing-power data for the charge-exchange reaction from PS199 |126| , |124| are shown in 
Figure |6|. Finally, in Figure |^ the prediction for the depolarization for elastic scattering at 
783 MeV/c is shown compared to the data from PS172 [ |150|| . Only a few depolarization data 
exist: one point at 679 MeV/c, three points at 783 MeV/c, and one point at 886 MeV/c. 
These data points are not included in the fit in view of the large error bars. This also means 
that no normalization error is taken into account. 

In Sect. VI a formalism was proposed to extract phase-shift parameters and inelasticities 
from the S matrix for antinucleon-nucleon scattering. It does not make much sense to present 
all these phase shifts, inelasticities and mixing parameters without a proper assessment of the 
uncertainties (statistical errors). This, however, requires a lot of work. Preliminary study 
shows that the phase-shift parameters for the 1 S'o and X P\ partial waves are not pinned 
down accurately at all above pi a b ~ 400 MeV/c. On the other hand, a large number of 
parameters appear to be very well determined by the existing data, such as the 3 -Po,i,2 an d 
3 Z} 2j 3 phase shifts and inelasticities, and the £x,2 mixing parameters. All results are available 
in numerical form from the authors. An extensive discussion of the remaining uncertainties 
in the predictions of the PWA and the phase-shift parameters will be presented elsewhere. 



In our 1991 preliminary PWA |2Pj we were able to determine the charged-pion-nucleon 
coupling constant / 2 from the data on the charge-exchange reaction pp — ► Tin, in which only 
isovector mesons can be exchanged. In this PWA we analyzed 884 scattering observables 
between 400 and 950 MeV/c, using 23 free parameters. We found / 2 = 0.0751(17), at the 
pion pole. The error is purely statistical. This value was in nice agreement with other deter- 
minations of / 2 from np [ 155[| and 7T p scattering data ||156|| . These results provided strong 
evidence for an approximate charge-independent pion-nucleon coupling constant, since they 
were consistent with the value for / 2 found in the Nijmegen pp PWA ||. 

Since this preliminary analysis, more high-quality analyzing-power data for the charge- 
exchange reaction have become available from PS199 ||124|| . We have repeated the deter- 
mination of / 2 , but this time from the complete 1993 Nijmegen database. The coupling 
constants of the neutral pion were kept at the value of fp = fl = 0.0745. Since / 2 is de- 
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termined by the data on elastic scattering pp — > pp, and / c by the data on charge-exchange 
scattering pp — > nn, one expects that the correlation between the two is small. Adding / c 2 
as the 31st free parameter, we now find 

f c = 0.0732(11) , 



at the pion pole. This result supersedes our previous value from Ref. ||29|| . Again, the error 
is of statistical origin only. In view of the enormous amount of work involved, it is very 
hard for us to make statements about possible systematic errors on this result. In Ref. |2S 
we did demonstrate that there were no systematic errors due to form-factor effects or due 
to p(770) exchange. In the Nijmegen pp PWA systematic errors could be more thoroughly 
investigated and there they were found to be small |K|. Although in our case the systematic 
errors are probably larger than for the pp case, it is very encouraging that the available 
charge-exchange data pin down the charged-pion coupling constant with a remarkably small 
statistical error and that the result is consistent with other determinations from 7r ± p ||156| 
and NN scattering ||155|,|60"1|. Very probably the new LEAR experiment PS206 ||153|| will 



further constrain the charged-pion-nucleon coupling constant. 



X. SUMMARY OF CONCLUSIONS 



To summarize, we have performed an energy- dependent partial-wave analysis of all 
antiproton-proton scattering data below 925 MeV/c antiproton momentum, published in 
a regular physics journal since 1968. This is the first time such an analysis has been at- 
tempted. We have set up the Nijmegen 1993 pp database by scrutinizing and passing sentence 
on all available pp scattering observables below 925 MeV/c. Serious problems were encoun- 
tered with a set of 540 elastic differential cross sections from LEAR and KEK. These data 
were rejected, although further study is required here. Of the remaining 3747 scattering 
observables 204 (5.4%) were rejected on the grounds of sound statistical criteria. We also 
rejected three normalization data. The final database contains 3543 scattering observables 
and 103 normalization data for a total of 3646 scattering data. Using 30 free parameters 
we obtain xLin = 3801.0 corresponding to Xmm/-^data = 1.043. This shows that the tail 
of the charge-conjugated Nijmegen potential is a realistic intermediate- and long-range pp 
force. Data on the charge-exchange reaction pp — > nn provide further evidence for a "low" 
and approximately charge-independent pion-nucleon coupling constant f^ Nn ~ 0.0745. The 
present results will serve as a starting-point for future investigations by our group of the 
antiproton-proton system, and should be of help in planning further experiments at LEAR 
and elsewhere. 
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APPENDIX: FROM S MATRIX TO PHASE PARAMETERS 

In his second paper on the subject of parametrizing the S matrix for nucleon-nucleon 



scattering above the pion-production threshold |94] , Bryan has presented an easy algorithm 
to extract the three parameters 5j-i,j, £ j, and the three different elements of the 

inelasticity matrix H 3 from the 2 x 2 S matrix, written as 

S J = exp(i5) exp(ie,ja x ) H J exp(ie,ja x ) exp(id) . (Al) 

For completeness we list the relevant expressions. The derivation can be found in the paper 
by Bryan. If the S matrix is presented numerically as 



s j ____ ( Ru exp(2z5 u ) iRi 2 exp(2iS 12 ) 
\iR 12 exp(2i5i 2 ) R22 exp(2i<5 2 2) 



(A2) 



then the phase shifts and <5j+i,j can be obtained from the following two equations 

R%r> sin 26 , . , 

tan2(0 a + 6 b ) = ^ , A3 

R11R22 + R12 cos ^ 

t&n(6 a - 6 b ) = R ™ ~ ^ n tan(^ a + 6 b ) . (A4) 

Here we defined the auxilary phases 

6 a = 5 U - 5j-i,j , (A5) 

e b = 5 22 - 8 J+1 ,j , (A6) 

5 = 5 n + 5 22 - 25 12 . (A7) 

The mixing parameter Ej can subsequently be calculated from 

tan 2, = ™^f*+J»-t) , (A8) 

R n cos 29 a + R 22 cos 29 b ' v ; 

where e = ej. Next the elements of the matrix H 3 are isolated. One finds 

2 cos 2e H u = R n cos 20 o (l + cos 2e) + R 22 cos 26 b (l- cos 2e) , (A9) 

2cos2efZ"22 = -Rucos20 o (l-cos2e)+i222Cos206(l + cos2e) , (A10) 

cos2sH 12 = R 12 sm(5-6 a -B b ) . (All) 



If the matrix H J is parametrized according to Klarsfeld [96], the three inelastic parameters 
Vj-i,Ji Vj+i,Ji an d U J can be obtained from 

7 lJ _ 1) j + 7 lJ+ ^j = r IiH J , (A12) 
V j_ lt jru +h j = det H J , (A13) 
tan 2toj = 2H 12 /{H n - H 22 ) . (A14) 
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FIGURES 



FIG. 1. Total and annihilation cross section as a function of momentum in antiproton-proton 
scattering. The otot data ctre from the LEAR, experiment PS172 |1 14|Jl"Q6| and. the <7ann 

data are 

from experiment PS173 [105, 10Cfl . The curve from the PWA for utot has Xmin = 88.4 for 75 points, 



and the curve for cr ann has x\ 



65.3 for 52 points. 



FIG. 2. Differential cross section for elastic scattering at 790 MeV/c. The data are from 



Eisenhandler et al. [|128| . The curve from the PWA has x mm = 101.5 for 95 points. 



FIG. 3. Elastic cross section at backward angle cos 9 = —0.994 as a function of momentum. 
The data are from Alston-Garnjost et al. [117]. The curve from the PWA has x min = 36.7 for 30 
points. 



FIG. 4. Differential cross section for charge-exchange scattering at 693 MeV/c. The data are 

Lmin 



from the LEAR experiment PS199 [126]. The curve from the PWA has Xmin = 39.3 for 33 points. 



FIG. 5. Analyzing power in elastic scattering at 544, 679, 783, and 886 MeV/c. The data 
are from the LEAR experiments PS172 [122] and PS198 [120]. The curves from the PWA have 



=37.5, 23.1, 30.1, and 38.5 for 27, 26, 29, and 34 points, respectively. 



FIG. 6. Analyzing power in charge-exchange scattering at 546, 656, 767, and 875 MeV/c. The 
data are from the LEAR experiment PS199 [126,124]. The curves from the PWA have Xmm =36.1, 



21.9, 25.5, and 20.0 for 23, 21, 22, and 23 points, respectively. 



FIG. 7. Depolarization in elastic scattering at 783 MeV/c. The data are from the LEAR 
experiment PS172 |15C|]. The curve from the PWA has x m in = 4-9 for 3 points. 



33 



TABLES 



TABLE I. Reference table of antiproton-proton scattering data below 925 MeV/c. 



Plab 


No. a 




Norm 


Pred. 








(MeV/c) 


type b 


v 2 . 

A mm 


error 


norm c 


Rejected d 


Ref. 


Comm. 



119.0- 










-923.0 


50 a cc 


29.2 


3-5% 


1.098 


176.8- 










-396.1 


5 <7 ann 


7.5 


4.4% 


0.943 


181.0 


46 dfiei 


48.7 


5% 


1.037 


183.0 


13 da ce 


8.3 


5% 


0.976 



194.8 

200.0- 
-588.2 

221.9- 
-413.2 

233.0 



19 dad 
48 <7 ann 
45 dtot 

54 d(7 e l 



57.8 

53.0 
87.6 



2.', 



70 



0.9% 

5% 



0.978 

0.971 
1.034 



239.2 


20 do- e i 


27.7 


4% 


1.089 


272.0 


65 d<7 e i 


55.9 


5% 


1.055 


276.0- 










-922.0 


21 cr cc 


32.0 


5-10% 


1.138 


276.9 


20 dad 


18.9 


4% 


1.042 


287.0 


54 da c \ 




5% 




287.0 


14 do- cc 


24.0 


5% 


1.201 


310.4 


20 d<7 e i 


30.1 


4% 


1.039 


340.9 


20 do- i 


32.2 


4% 


1.044 


348.7 


38 do- e i 


42.4 


4% 


0.993 


355.0- 










-923.0 


36 otot 




1.5% 




353.3 


119 do- c i 


116.4 


5% 


1.037 


369.1 


19 do-ei 


15.7 


4% 


1.020 


374.0 


39 da c \ 


24.3 


5% 


1.067 


388.0- 










-598.6 


29 ot t 


35.4 


0.7% 


0.973 


392.4 


19 do-ei 




5% 




392.4 


15 do- cc 




5% 




404.3 


40 do-ei 


38.6 


4% 


0.994 


406.0- 










-922.0 


30 do- e i 


36.7 


4% 


0.896 



176.8 
>0.925, #=6 
0.940, -0.170, 
-0.574, -0.966 
all 



221.9, 229.6 
>0.938, #=6; 
0.764 



0.967 



all 

0.985 

-0.850 



all 
0.366 
0.550 



all 
all 

-0.575, -0.925 



[100| 



L00,105 



HI 



[104| 



mm 



mm 



HioTl , |To2 



HI 

pi 



[110] 



E3 



EH 



[116] 



k, m 



J, o 



f, i, o 



J 
o 

J 

m 

o 

j, 1, o 

o 
o 
i, o 

e, m 
j, o 
i, o 
o 



1 
f 

i, o 
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TABLE I. (Continued.) 



Plab 


No. a 




Norm Pred. 








(MeV/c) 


type b 




error norm c 


Rejected d 


Ref. 


Comm. 



406.0 
411.2 
413.4 
424.5 
428.0 
435.8 
439.0 
439.0 
439.9 
440.8 
444.1 
446.0 
447.1 
458.3 
467.5 
467.8 
469.2 
479.3 
480.0 
481.2 
490.1 
490.1 
490.6 
492.7 
497.0 
498.7 
503.8 
504.7 
505.0 
505.0 
508.0 
508.9 
516.0 
523.0 
524.8 
525.9 



119 d(Tel 

38 ddei 
7 dad 
7 do-ei 

10 docc 
7 do-ei 

27 do-ei 

24 A Viel 

39 do- e i 
38 do- e i 

38 do- e i 

119 dCT e i 

7 do-ei 

8 do-ei 

39 da e i 
39 da e i 
8 do-gi 

119 do ol 
10 dOce 
8 do-ei 
37 do- e i 

15 dOce 

39 do- e i 
8 do-ei 

14 Ay >el 

37 do- e i 
8 do-ei 
39 do- e i 
54 do-gi 

14 dOce 

119 do cl 
39 da e i 
8 do-gi 

15 Ay je \ 

36 do- e i 
39 do- c i 



99.8 
37.9 
5.6 

12.3 
1.6 

38.8 
42.0 
61.4 
35.0 
115.3 

6.9 

2.0 
32.7 
24.0 

8.2 
109.3 
14.3 

7.2 

12.5 
47.5 
4.2 
7.3 
27.7 
13.3 
14.3 

30.1 
105.3 

28.7 
5.5 
8.3 

32.2 

45.0 



o/o 
5% 
5% 
5% 
20% 
5% 
10% 
5% 
5% 
5% 
4% 
5% 
5% 
5% 
4% 
5% 
5% 
5% 

oo 
5% 
5% 
5% 
5% 
5% 

oo 
4% 
5% 
5% 
5% 
5% 
5% 
5% 
5% 

oo 
4% 
5% 



1.029 
1.025 
1.054 

1.221 
1.017 

1.068 
1.031 
1.035 
0.972 
1.021 
1.050 
0.994 
1.039 
1.056 
1.013 
1.003 
1.154 
1.048 

1.068 
0.983 
1.014 
0.718 
1.004 
1.047 
1.021 

1.034 
1.019 
1.025 
1.018 
0.786 
1.020 
1.053 



0.990, 0.750, 0.578 
-0.925 
0.992 
all 

0.992 

all 
0.851 



0.175, -0.825, -0.875 

0.992 
0.996 
-0.925 

0.996 
0.919, 0.873, 0.697 

0.996 

all 
-0.193 

0.996 



0.996 

all 
0.574 
0.663, 0.530 

0.996 



[112| 



[113 1 



[118| 



[119| 



[118| 



[120| 



113] 



113 1 



110] 



112] 



[1181 



[118| 



110 



113| 



[118| 



121] 



118 1 



116 



113 



118 1 



122,123] 



110 



118 1 



113 



[|101| , |102 | 



103 



[112 1 



113 



118 



122,123] 



110 



113 



J) o 
i, o 
j, o 
e, j, o 

j, o 
1 
o 
o 
i, o 
o 
o 
o 
o 
1, o 
o 

j, o 
j, o 

g 

j, o 

1 



J) o 

h 

i, o 
j, o 
o 

j, 1, o 
j, o 



j, o 

h 
i, o 

o 
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TABLE I. (Continued.) 



Plab 


No. a 




Norm Pred. 








(MeV/c) 


type b 




error norm c 


Rejected d 


Ref. 


Comm. 



528.2 
533.6 
537.0 
540.6 
543.2 
544.0 
544.0 
546.0 
549.4 
550.0 



8 dcr e i 

119 dCTel 

10 dcr cc 
8 dcr el 
39 d<T e i 
33 ddei 
30 A y>el 

23 Ay^QQ 

10 dace 
67 dcr c i 



3.2 
133.7 
19.4 
10.9 
38.4 

37.5 
36.1 
7.1 
76.0 



o/o 

5% 

oo 

5% 

5% 

10% 

5% 

4% 

20% 



1.005 
1.029 
1.199 
1.015 
1.065 

1.046 
0.959 
1.258 
1.006 



553 1 


34 drr i 

•) i LLC q| 


38 6 


4% 


981 


553.4 


8 drr r 


2.4 


5% 


1.017 


556.9 


119 do- i 


125.4 


5% 


1.025 


558.5 


39 da i 


45.4 


5% 


1.040 


565.5 


8 drr i 


5.7 


5% 


1.006 


568 4 


37 Hrr i 


34 3 


5% 


1 040 


577.2 


36 do"ei 


36.0 


4% 


0.983 


578.1 


9 dflel 


6.2 


5% 


1.014 


578.3 


119 da el 


132.3 


5% 


1.047 


584.0 


10 do- co 


15.7 


oo 


1.112 


590.0 


39 ddei 




5% 




590.0 


15 dflce 


32.8 


5% 


1.092 


591.2 


9 ddei 


6.5 


5% 


1.029 


591.2 


39 ddei 




5% 




591.2 


15 dflce 


18.0 


5% 


1.058 


596.5 


38 da e \ 


49.5 


5% 


1.075 


599.2 


33 ddei 


15.8 


4% 


0.997 


604.0 


9 da- e i 


7.6 


5% 


0.987 


615.0 


38 da e \ 


48.1 


5% 


1.046 


617.0 


9 dcr e i 


6.7 


5% 


0.953 


630.0 


10 dflce 


9.3 


oo 


1.073 


630.9 


9 d(J e l 


4.6 


5% 


0.997 


639.6 


38 ddei 


17.1 


5% 


0.995 


644.7 


9 do- c i 


7.8 


5% 


0.996 


656.0 


21 Ay Ce 


21.9 


4% 


0.963 


658.1 


38 d<T e i 


37.4 


5% 


0.972 



0.996 



0.996 
-0.975 
all 

>0.883, #=3 



>0.995, #=3; 
0.910, 0.883, 0.869 

0.996, 0.972 
0.908 

0.996 
-0.675, -0.825 

0.999 



all 

0.996, -0.574 
0.999 

all 
-0.358 



0.998 
-0.575 
0.998 

0.998 
-0.175 
0.998 

0.225, -0.675, -0.975 



[118| 



[112 1 



[118| 



[113| 



[120| 



pg 



124 1 



119 



|25J 



110| 



[1181 



[H2| 



113 1 



[118| 



113 1 



HQ I 



112| 



121| 



H1011JL02 I 

H 



118 1 



[115| 



116 



113 



110 



118 1 



113 



118 1 



121 



118 1 



113 
118 




124 


,126 


113 





J) o 

j, o 

g 

j, o 
o 

1 

o 



J 

i, o 
j, o 
j, o 
o 

, o 
, o 
, o 
, o 

J) o 



j, 1, o 
j, o 

1 

o 
i, o 
j, o 

o 

j, o 

g 

j, o 
o 

j, o 
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TABLE I. (Continued.) 



Plab 


No. a 




Norm 


Pred. 








(MeV/c) 


type b 


X 2 - 

Amin 


error 


norm c 


Rejected* 1 


Ref. 


Comm. 



658.6 


9 dov.1 

'■' A v el 


8.9 


5% 


1.004 


0.998 


670.0 


10 der„„ 


5.5 


oc 


1.165 




671.5 


9 d(Tnl 


3.3 


5% 


0.987 


0.998 


679.0 


26 dcr e i 




OC' 




all 


679.0 


27 -^y^el 


23.1 


OC' 


0.846 


0.540 


679.0 




1.4 








679.1 


4 A, „i 


6.3 


5% 


0.983 




680.1 


38 dcr c i 


40.2 


5% 


1.003 




686.1 


9 dcToi 


3.9 


5% 


0.986 


0.998 


689.0 


39 dcr c i 




5% 




all 


689.0 


16 dcr„„ 

1 ' 1 ' cc 


26.5 


5% 


1.010 


-0.139 


690.0 


89 dcr e i 


103.5 


4% 


0.991 


0.370 


693.0 


34 dcr„„ 

W VAV^ {^g 


39.3 


10% 


1.069 


-0.075 


696.1 


21 davi 


16.4 


4% 


1.026 




696.1 


16 d<7™ 

V/ VAU (^g 


21.0 


4% 


1.050 




697.0 


24 d(j P i 




10% 




all 


697.0 


33 g} 


20.3 


5% 


1.022 


0.629 


698.0 


10 dcT™ 


7.1 


oc 


1.237 




700.0 


4 A, pi 


2.1 


5% 


0.991 




701.1 


9 dfjoi 

t/ g^ 


3.5 


5% 


1.000 


0.998 


715.3 


9 dfjoi 


10.6 


5% 


1.002 


0.998 


728.0 


10 d<r™ 


2.9 


OC' 


1.105 




757.0 


72 dcj P i 


95.7 


5% 


1.055 


>0.996, # 


767.0 


22 A, r P 
j A y,ce 


25.5 


4% 


1.113 




780.5 


39 d<7 e i 




5% 




all 


780.5 


15 dcr^„ 

v A ' ce 


14.0 


5% 


0.974 




783.0 


30 d<7 e i 




OC' 




all 


783.0 


30 A, el 


30.1 


4.5% 


0.944 


-0.300 


783.0 


3 Dyy 


4.9 








790.0 


95 dcr e i 


101.5 


4% 


1.034 




860.0 


95 dcr e i 


70.5 


4% 


1.045 


0.510 


875.0 


23 Ay Ce 


20.0 


4% 


0.995 




886.0 


34 do-el 




CO 




all 


886.0 


34 Ay jS i 


38.5 


4.5% 


0.992 




886.0 


1 Z) w 


1.2 








910.0 


19 dcr c i 




CO 




all 


910.0 


21 ^4y,el 


14.7 


5% 


0.989 





1 118 



1 121 



m 



|122,123| 



1 150 



127 



1 113 



Big 



1 115 



| 116 



128 



1 126 



[129| 



[130| 



[121 



HI 

[ 118 | 
1 118 1 

m 

1 125 1 
[124| 

\m\ 

[ 116 1 

a 
muni 

m 

[12|l 
[128| 

[124| 

m\ 

[ I22]|123 | 

in 

[1J2| 

PHI 



J) o 

g 
j, o 

1 

h 

p, q 

o 
o 

j, o 

1 



1 

o 

g 

o 

j, o 
j, o 



1 
1 

p, q 



p, q 

f, e 
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a The number includes all published data, except those given as 0.0±0.0 (see Comment 
i), and those having pi ab > 925 MeV/c (see Comment m). 

b The subscripts "el" and "ce" denote observables in the elastic pp — > pp and charge- 
exchange pp — > nn reactions, respectively, "do - " denotes a differential cross section 
da/dQ, A y a polarization-type datum (asymmetry or analyzing-power), and D yy a 
depolarization datum. <7 to t stands for total cross section, a ann for (total) annihilation 
cross section, and a cc for integrated charge-exchange cross section. 

c Normalization, predicted by the analysis, with which the experimental values should 
be multiplied before comparison with the theoretical values. 

d Tabulated is pi a b in MeV/c or cos 9. The notation ">0.925, #=6" e.g. means that the 
6 points with cos 9 >0.925 are rejected. 

e Group rejected due to improbable low Xmm- 

f Group rejected due to improbable high Xmm- 

g Floated normalization. Data are relative only. 

h Normalization floated by us, since the norm contributes much more than 9 to xLm- 
i Data points given as 0.0±0.0 not included. 

j Coulomb- nuclear-interference measurement. Data points in the extreme forward an- 
gular region are rejected when they contain multiple-scattering effects. 

k Data points at low momenta rejected (see text). 

1 Problematic differential cross section. Not included in the database. See the text, 
Sect. VIIIB and Tables II and III. 

m Part of a group of data with points having p lab > 925 MeV/c. 

n Elastic differential cross sections as a function of momentum taken at backward angle 
cos 9 = -0.994. 

o Normalization error assumed by us, since no clear number is stated in the reference, 

p Depolarization data. Not included in the fit, in view of large error bars, 

q Normalization error taken to be zero, in view of large error bars of these data, 

r Data points taken at the same angles averaged. 
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TABLE II. Corrections applied to elastic differential cross sections from Eisenhandler et al, 
from PS172, PS173, and PS198 from LEAR, and from KEK. The number of data includes all 
published points. The column labeled "Syst. error" gives the approximate point-to-point error 
that has to be added quadratically to the statistical error to reach Xmm ~ ^df- The column 
labeled "Xunc" gives the result of the Legendre fit without any corrections, except that PS173 
points at forward angles contaminated by multiple-scattering effects are removed. For the values 
of Xmm after the corrections, see Table P4 



Plab 


No. 


Syst. 


(MeV/c) 


Group data Rejected points (cos 9) 


error Xunc ^ e ^- Comm. 



181.0 


PS173 


46 


287.0 


PS173 


54 


392.4 


KEK 


19 


439.0 


PS198 


27 


490.1 


KEK 


18 






19 


505.0 


PS173 


54 


544 

Oil. U 


PS1 98 




590.0 


PS173 


39 


591.2 


KEK 


19 






20 


679.0 


PS172 


26 


689.0 


KEK 


15 






24 


690.0 


Eisenh. 


89 


697.0 


PS198 


24 


780.5 


KEK 


14 






25 


783.0 


PS172 


30 


790.0 


Eisenh. 


95 


860.0 


Eisenh. 


95 


886.0 


PS172 


34 



> 0.92, #=6 

> 0.95, #=9; 0.345, 0.199, 0.101, 0.051, 



-0.345 



> 0.98, #=3; 0.96, 0.67 

> 0.99, #=2 
0.825 

0.275, -0.075 
0.50, -0.50 



0.37 



0.58, -0.58 



0.66, 0.54, 0.46, -0.18, -0.66 



5% 
2% 
4% 
4% 



41.0 
160.7 
19.2 
28.0 
202.0 
15.0 
59.5 
23.1 
39.0 
276.2 
36.0 
82.7 
90.5 
24.2 
101.9 
174.0 
30.5 
16.4 
31.9 
88.2 
69.4 
77.6 



101.102| 



lQgJTogl 



11511 



he] 



iQj]J102 | 



10l|Jl02 | 



115] 



123] 



115] 



125 



130] 



115] 



123H 



125 



125 



122] 



a Data points in the extreme forward angular region are rejected because they contain 
multiple-scattering effects. 

b Data in the "one-prong" region. For a discussion, see the text. 

c Data in the "two-prong" region. For a discussion, see the text. 

d The most forward and the most backward data point are rejected (F. Bradamante, 
private communication). 
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TABLE III. Legendre-polynomial fits to elastic differential cross sections from Eisenhandler et 
al, PS172, PS173, and PS198 from LEAR, and from KEK. The corrections summarized in Table |Tj 



are taken into account. The number of data is the number of published data minus the rejected 
points. For the KEK groups at 490, 591, 689, and 781 MeV/c the first line gives the results for the 
"one-prong" region, and the second line the results for the "two-prong" region. 



Plab 

(MeV/c) 


Group 


No. 
data 


a 


oi/cio 


a 2 /a 


a 3 /a 


a 4 /a 


a 5 /a 


a 6 /a 


2 

Xmin 


181.0 


PS173 


40 


7.02(25) 


0.92(06) 


0.46(07) 










41.0 


287.0 


PS173 


40 


5.10(07) 


1.58(03) 


0.82(03) 


0.15(02) 








36.0 


392.4 


KEK 


19 


5.49(06) 


1.92(04) 


1.35(05) 


0.35(03) 








15.6 


439.0 


PS198 


27 


4.09(08) 


1.93(06) 


1.63(07) 


0.67(08) 


0.39(09) 


0.19(07) 


0.09(04) 


19.6 


490.1 


KEK 


18 


5.07(10) 


1.95(05) 


1.55(07) 


0.50(05) 


0.10(03) 






15.4 






19 


5.61(45) 


1.91(21) 


1.89(28) 


0.72(13) 


0.38(10) 






15.0 


505.0 


PS173 


49 


4.18(05) 


1.90(04) 


1.76(04) 


0.96(05) 


0.49(06) 


0.11(05) 


0.06(04) 


36.1 


544.0 


PS198 


33 


3.79(05) 


2.08(05) 


2.16(06) 


1.16(06) 


0.57(06) 


0.12(04) 


0.04(02) 


23.1 


590.0 


PS173 


37 


3.73(06) 


2.11(05) 


2.09(05) 


1.25(06) 


0.66(06) 


0.26(05) 


0.10(03) 


32.2 


591.2 


KEK 


18 


4.42(08) 


2.12(05) 


1.95(06) 


0.83(04) 


0.16(03) 






13.8 






18 


5.24(33) 


2.07(18) 


2.08(22) 


0.99(11) 


0.38(07) 






15.2 


679.0 


PS172 


24 


1.96(02) 


2.59(05) 


1.88(02) 


1.39(04) 








18.8 


689.0 


KEK 


15 


3 52('22 N > 


2 27f1 51 


2 35f1 


1 46f4 1 1 


631441 


1 7(07) 




q 8 






24 


3.21(15) 


2.16(14) 


2.14(17) 


1.20(09) 


0.38(06) 






19.0 


690.0 


Eisenh. 


88 


4.28(07) 


2.28(06) 


2.45(08) 


1.67(07) 


0.76(06) 


0.23(04) 


0.03(02) 


81.6 


697.0 


PS198 


24 


3.54(09) 


2.22(08) 


2.30(09) 


1.43(08) 


0.51(05) 


0.10(02) 




21.5 


780.5 


KEK 


14 


3.14(23) 


2.37(17) 


2.72(23) 


2.10(16) 


1.07(16) 


0.29(07) 




9.1 






25 


2.39(10) 


2.07(11) 


2.39(16) 


1.36(08) 


0.66(07) 






16.4 


783.0 


PS172 


28 


2.55(07) 


2.73(17) 


2.45(11) 


2.26(20) 


0.67(05) 


0.38(07) 




29.8 


790.0 


Eisenh. 


95 


3.87(08) 


2.36(07) 


2.69(09) 


2.08(09) 


1.09(07) 


0.42(05) 


0.12(03) 


88.2 


860.0 


Eisenh. 


95 


3.67(06) 


2.41(06) 


2.83(07) 


2.33(07) 


1.34(05) 


0.56(04) 


0.17(02) 


69.4 


886.0 


PS172 


29 


2.44(04) 


2.66(08) 


2.64(06) 


2.46(09) 


0.95(03) 


0.47(04) 




22.2 
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TABLE IV. P-matrix parameters of the different partial waves. Vo and V\ are the real parts 
of the short-range spherical-well potential, for isospin 1 = and 1 = 1 respectively. W is the 
isospin- independent imaginary part. The mixing angles that take care of the short-range 1 = 
tensor force are: 6>i = 55.7° ± 1.3° for the 3 Ci state, 9 2 = 45.8° ± 1.1° for the 3 C 2 state, and 6 2 = 
10.7° ± 4.7° for the 3 C 3 state. The quoted errors are defined as the change in each parameter that 
gives a maximal rise in Xmm °f 1 when the remaining parameters are refitted. 



partial wave 


V (MeV) 


Vi (MeV) 


W (MeV) 




n 

u 


n 

u 


— yy^u j 




-151(6) 


-17(3) 


-100(3) 


'Pi 








-90(7) 




-132(9) 


178(19) 


-156(9) 


3 P1 


155(12) 


-66(3) 


-97(4) 


3 P2 


-136(4) 


-69(3) 


-142(3) 


l D 2 








-105(12) 


3 D, 


-215(7) 


33(16) 


-106(5) 


3 D 2 


-38(12) 


-198(4) 


-110(4) 


3 D 3 


-152(6) 


-102(5) 


-163(4) 


3 F2 


-101(20) 


-250(14) 


-179(7) 
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TABLE V. Partial elastic and charge-exchange cross sections in mb. 







pp 


->■ pp 






pp 


— ► nn 




Plah (MeV/c) 


200 


400 


600 


800 


200 


400 
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